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An integral of the Wigner function of a wavefunction \ip >, over some region 
S in classical phase space is identified as a (quasi) probability measure (QPM) of 
S, and it can be expressed by the \ijj > average of an operator referred to as the 
region operator (RO). Transformation theory is developed which provides the RO 
for various phase space regions such as point, line, segment, disk and rectangle, and 
where all those ROs are shown to be interconnected by completely positive trace 
increasing maps. The latter are realized by means of unitary operators in Fock 
space extended by 2D vector spaces, physically identified with finite dimensional 
systems. Bounds on QPMs for regions obtained by tiling with discs and rectangles 
are obtained by means of majorization theory. 



Key words: Wigner function, region operators, Quantum Mechanics, quantum tomog- 
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1. Introduction 

The Wigner quasidistribution function on phase space W(p, q) [1], provides an impor- 
tant tool in the description of quantum systems from both theoretical and experimental 
point of view. Of particular interest is for example, its application to the field of quantum 
tomography [2] , [3] . Although the Wigner function shares several properties with classi- 
cal distribution functions on phase space it is not positive-definite. As a consequence, 
the value of its integral, (to be called alternatively quasiprobability integral(qpi), mass 
(qpm), or volume (qpv)), over subregions of phase space, lies in general outside the inter- 
val [0, 1], and may admit large positive or negative values. The upper and lower bounds 
of such integrals define tests e.g on the accuracy of experimental determinations of the 
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Wigner function in the sense that they give the degree to which the value of the integral 
lies outside [0, 1], and thus their determination is of particular importance. In [4], [5], 
[6] and [7] it has been shown that such bounds, for given regions or contours of phase 
space for ID quantum systems, are determined by the maximal and minimal eigenvalues 
of corresponding Hcrmitian operators K$, called region operators (RO), whose spectral 
analysis has been shown to have a rich algebraic structure. Consequently the study of 
the properties of such operators defined by 

K s = J^ X s(a)W(a)da , (1) 

where S is a subregion in phase space r, a <E T, xs( a ) is the characteristic function 
corresponding to S and W^a) are the Wigner operators, is of great importance. 

This article presents a detail analysis of transformation theory for RO, and in par- 
ticularly investigates the cases of the point, line, rectangle and disk operators. The 
development of the theory is mainly based on the notion of completely positive trace 
increasing map (CPTI) (c.f [8]), which provides the interconnection between various re- 
gion operators. The application of such maps is also shown to describe tiling processes in 
phase space. Moreover the eigenvalues of the obtained ROs from step to step during the 
tiling are shown to respect majorization relations[10]. This allows the determination of 
upper and lower bounds on qpi over regions in phase space obtained by a tiling process. 

An outline of the work is as follows: section two focuses on line and rectangle ROs. 
The spectral analysis of straight line segment ROs is presented and further it is shown that 
straight line ROs are projection operators. On the other hand the normally ordered form 
of rectangle ROs is obtained, and moreover it is shown that all rectangle ROs obtained 
from an initial one by rigidly shifting its position in the phase plane, are isospectral and 
thus have the same quasiprobability volume. 

In section three an operational construction of quasiprobability measures of canonical 
polygon ROs is developed by introducing appropriate completely positive trace increasing 
(CPTI) maps. An explicit example is then given in which a general canonical hexagon 
RO is constructed from the corresponding RO of an isosceles triangle, by means of CPTI 
maps. 

Finally in section four, appropriate CPTI maps are used to realize a tiling process in 
phase space acting on rectangle and disk ROs used as building blocks. The step matrix 
which connects the eigenvalue- vectors of ROs at each step of the tiling process is shown 
to be simply determined by doubly stochastic matrices. This further allows to determine 
simple ordering relations between the eigenvalues of successive ROs by means of the 
majorization theory. This fact provides intervals of extremal values for the eigenvalues 

1. e the qpm's, of the sequence of region operators, in a tiling process. Section five briefly 
summarizes results and outlines some prospects of the theory of region operators. 

2. Point, Line and Rectangle Region Operators 

Proposition 1. Region operators having support on a straight line segment of length 
L in angle 9 with the q—axis read, K e L = sln ^Qc°se+PsinP L ^ an d satisfy the eigenvalue 
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problem K e L \^ e ± ) = ± 2sin ^ eL) \V e ±), with eigenvectors \^ e ± ) = \q e ) ± \-q ), where {\q g ), 
qg G R}, is a complete set of vectors determined in the number state basis by means of 
the Hermite polynomials as (qe\n) = 2^/2^1)1/2 cx P(~ge/2)-Hn(ge) exp(— wfl). 

Proof: Let K be the region operator with support on the point at the origin ("point 
operator", see fig. la), identified with the parity operator II, since 

Kq = J R2 S(q)S(p)D(q,p)HD(q,p)^dqdp = II. Let us introduce a symmetric broadening 
of the point at the origin along the q— axis of extension L, by means of the positive 
map K£ = e^(Ko) = e£(II), constructed from the displacement operator D(q,p) = 
expi(pQ — qP) — exp(ao^ — a*a), with a = -^(q + ip) and a = -^{Q + iP), a) = 
-j=(Q-iP), to get 

— — 

K P L = f \ D{q, 0)110(9, 0)Uq = f \ D{2q, 0)dgII = S J^1 U . (2 ) 

In this way we obtain the region operator K£, with support on the straight line segment 
of length L, symmetric with respect to the origin and lying along the position axis, (see 
figure lb). Similarly for the segment along the momentum axis the operator is = 
£ Q(K ) = SHiiipn. Since U\q) =\-q), we get K^f) = ±^Md\^Q) 7 where = 
\q) ± I — q), q € R. In the same manner for the region operator with support a on length 
L line segment centered at the origin and extended along the p— axis, we get K^\^f^.) = 
_l_ 2sin(pL) |v|/P^ where = \p) ± | —p),p € R. The generalized eigenvectors, have 
respectively the following position and momentum representation: ip±(q ) = (q \^±) = 
5(q +q)±6(q - q), and j>%') = (p'|*«) = ^ cos(«p'), i> Q (p) = {p'\^ Q ) = 
7=sin((7p'). 

Similarly for the position and momentum representations of the generalized eigenvec- 
tors we correspondingly obtain the functions ) = (Q |*+) = ^= cos(pq ) ,ip^(p ) 

= ^ sin ( w ')' and ^±( p ') = (p'I*±) =<5 (P +p) ±(5 (p' -P)- 
If we now rotate by angle 9 the line segment support of region operator K®, by the 
rotation operator e t6N , where AT the number operator, we end up with the operator 



= J0N K Q -UN = sin[(Qcos6> + Psing)L] n 
Ai L Qcos6 + Psm9 



As special cases we obtain the region operators for the line segment along the p— axis 
K 8=o = K Q^ and along the q _ axig K e=*/2 = r p u 

Remarks: 1) As the trace of a region operator equals the area of the support of the 
operator itself, the two positive maps e^' P : Kq — > K® P , introduced previously to get 
the straight line segment operator from a point operator are trace increasing maps i.e 
Tr{Kl) = Tr(e Q L (K )) = LTr(K ). 

2) The region operator K L = sin (Q<> L ) n, is written in terms of the rotated operator 
Qg = QcosO + Ps'm8. Two copies of the latter with their angles of rotation differing by 
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tt/2, 

Qe — Qcosd + PsinO, Qe+n/2 — Qsin# - Pcosd, (4) 

are canonical i.e [Qe, Qe+-n/2] = il- 

3) Some additional properties of the family of vectors {\qe), qe € R}, are: it is a 
complete set i.e dgelgeX^el = 1, its generalized vectors are orthogonal i.e (qe\q g ) = 
6(q$ — q' g ), and can be constructed from the zero number state as 



i r i i 



t2 



|0>, (5) 



and so the parity operator flips the sign of its argument i.e H\qg) = \ — qe)- These same 
states are also specified from their overlap with a coherent state i.e 

1 1 i 

(le\a) = ^T7i ex P [i{Qe+ir/2l<le\ exp[--(q - {Qe)) 2 ] exp[--(Q e )(Q e+ff / 2 )], (6) 

where {Qe) — {&\Qe\ot) — ^[aexp(— id) + a* exp(iB)]. 

4) Region operators with support on straight line segments with proportional lengths, 
in any direction, are commuting i.e K 6 XL K e L = K e L K^ L , for A > 1, and < 6 < 2n. 

PROPOSITION 2. Region operators having support on a straight line are projection 
operators. 

Proof: Let the region operator Kq = II. Its smearing along the whole momentum 
p— axis results into the following region operator 



/OO 1 
d P D(0,p)ILD(0, P y = -= : e-Q 2 := \q = 0){q = 0|. 
-oo V 71 " 



(7) 



If we next transform this operator to have support along an axis parallel to p— axis and 
crossing the position axis in the point q £ R then we will obtain 

K Q=q = D(q,0)K Q=0 D(q,0y = -L : ^ q -^ := \q){q\, (8) 

namely the projection operator in position states. Similarly by smearing the point oper- 
ator along the position q— axis we obtain 

/OO i 
dqD(q, 0)UD{q, 0)t = — : e - p2 := \p = 0) (p = 0| . (9) 
-oo V 71 " 

This operator meets the p— axis at its zero point, so by displacing it arbitrarily by p G R 
we will obtain 

K P=p = D(0,p)K P=o D(0,p)l = 4= : ^ (p - Pf ■= |p><p|, (10) 

V 71 " 
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namely the projection operator in momentum states. In the last four equations above 
the expression of respective region operators has been given also in normal ordered form 
i.e when the canonical creation operator is placed in the left of the annihilation operator. 

So far operators Kq =q = \q){q\, and Kp =p — \p)(p\, which are projective and have 
support on any straight line parallel to the p— axis crossing the q— axis at the point q, and 
respectively on any straight line parallel to the q— axis crossing the p— axis at the point p 
have been constructed. In order to rotate these support axes by any desired angle 8 with 
respect to e.g the q— axis, we should perform an additional unitary rotation with e l6N , 
where N is the number operator, the generator of rotations around the origin. Namely, we 
should consider the transformation K Q=q (K P=p ) -» e WN K Q=q e~ i0N {e i6N K P=p e- i6N ). 
The latter will not change the projective character of Kq =q = \q)(q\, and Kp =p — \p){p\, 
operators^ (see also [9]) 

Remarks: 1) The region operator having support on the whole plane is the unit 
operator. In order to obtain the region operator on the plane we should "add" the region 
operators on all lines parallel to any given line. Let e.g Kq =q = \q)(q\, and Kp =p = \p}{p\, 
the two projective operators that have support on some straight line parallel to p— axis 
crossing the q— axis at the point q, and respectively on some straight line parallel to the 
(j— axis crossing the p— axis at the point p. Their integrals J R KQ =q dq = J R \q)(q\dq = 1, 
and J R Kp =p dp = J R \p)(p\dp = 1, give the region operator supported on the plane, 
which actually is the unit operator due to the completeness of position or respectively 
momentum states. The same result would have been obtained if we would had started 
with any line in angle 9 with the p, q axes. 

2) Various region operators can easily be introduced by combining the operators 
constructed above e.g a region operator with support on a bundle of n parallel lines 
crossing the q— axis at the points P n = {q\, q^.--, q n } is defined as K Pn = \q\){q\\ + 
|f?2) {'22 1 + ••• + |<7ra)(9ra|- The latter can be obtained acting on |<7i)(gi|, the operator with 
support on the first line of the bundle, with the trace increasing map s n , i.e Kp n = 
e n (\qi)(qi\) — Y^i=i i which has an operator sum representation in terms of 

the unitary operators V% = exp(i(qi — qi)P), i — 1,2, ...,n. 

Let us now show a more special result related to motions that can be done on rectan- 
gle region operators that leave their spectrum, namely their associated quasiprobability 
volumes, invariant. 

Proposition 3. If we denote by K$ = K(xo, k n ; A, B) the region operator with 
support on a rectangle domain S in phase space with lower left corner at the point 
(xo, ho) and sides of lengths A and B along the axes of position and momentum variables 
respectively, then i) the normal order form of this operator is given by equation (14) ii) 
all operators K(xo, fco! A, B), K(xq + s, ko + 1; A, B) for s, t G R are isospectral, namely 
they are region operators obtained from the initial rectangle operator K(xo,ko; A, B) by 
rigidly shifting its position by s,t, and have all eigenvalues equal, which then implies that 
they have the same quasiprobability volumes. 

Proof : i) Let us recall that in general an operator Ks associated to a phase space 
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region S is given in terms of the Wigner operator 



W(a) = ±D{a){-l) N D(a)1 = ±e 



\ rw , J _ J_ p iTr(a t -a*)(a-a) _ _J_ . e -2(a< -a")(a-a) . ^-q) 



2tt 



The last equation is in normally ordered form and it is easily obtained by means of the 
formula e A ( at ~" )( a ~ a ) =: e( e * _1 )( at_a )( a ~ Q ) : valid for complex parameter A. Then by 
means of the (over) complete basis of coherent state vectors \z = q + ip) = \qp), z e C, 
the region operator becomes 



-I I 

2?r J zee J a 



(2tt)2 



zeC Jaesupp(xs) 



. e -2(at- Q *)(a-a) d 2 a . 

zeC Jaesupp(xs) 

e -2(z*- a *)(z-cx)\ z ^ z \ d 2 a d 2 z . 



(12) 



At this point we employ the technique of integration within an ordered product 
(IWOP)[ll], which applies in our case since the above operator integrand is in normally 
ordered form. In effect splitting the integration and applying IWOP yields, after using 
the integral representation of error function erf(x) = J* e~* dt , the normal ordered 
expression of the region operator 



K s = 7A2 : / ( / 
(2tiT J(q, P )eR 2 Jx 



Xq+A 



-{q-xY 



dx){ 



1 

16tt 



J ko 



k +B 



-(p-ky 



dk) \qp) {qp\dqdpdxdk : (13) 



er/(--^=(a t + a) + x + A) - er/(--^=(a t + a) + x ) 
er/(--^=(a t - a) + ko + B) - er/(--^=(a t - a) + k ) 



: . (14) 



ii) If we rewrite the rectangle operator Ks = K(xq, ko; A, B) in terms of position 
and momentum operators i.e Q = -^{^ + a,), P — "^( a ^ ~ a )> usc the displacement 
operator property £)(»^)QD(2^)t =Q-s, D(^)PB(^f ) f - P-t, and the series 
expansion of the error function, we obtain that 



D(^ 7 ^)K(x ,k :A,B)D(- 



it 



V2 



V2 



) ! = K(x + s,k +t;A,B). 



(15) 



Interpreting these two formulas we say that the family of region operators resulting 
from rigidly shifting the vertices of the initial rectangle support i.e S(xq, ko; A, B) — > 
S (xo + s, ko + t; A, B) , associates itself with region operators that are unitarily equivalent 
to the region operator of the original rectangle K(xo, ko; A, B). This implies that all these 
operators have the same spectrum and in turn means that their respective eigenvalues 
i.e the quasiprobability volumes are equal ■ 
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3. Operational Construction of QPMs: The Case of Canonical Polygons 

Let us introduce rotation operators in terms of the number operator N, to be the 
operators R(^f) = cxp(^-jiV), which generate rotations by angles j = 0, 1, M — 
1. Let us introduce an M— sided canonical polygon with radial parameter a, defined as 
the distance between the center of the polygon and the midpoint of a side. Due to its 
geometric invariancc under rotation by angles -jff-, j = 0,1,.. .,M — 1, the associated 
polygon region operator denoted by X[ 0) m] is obtained as the sum of copies of the region 
operator X^ a ^ of an isosceles triangle, rotated by all angles j = 0, 1, M — 1, [7], 
i.e 

*(*[«,*]) - *[.,*] = E^)*!-,*]*^)*- (16) 

3=0 

In the above equation the completely positive trace increasing map (CPTI) e, has been 
introduced by means of its M generators, given by the rotation operators. Since the trace 
of a region operator equals the area of its associated region, the CPTI map e that produces 
an M— polygon operator from its generating isosceles triangle region operator, should be 
a trace (area) increasing one. Indeed we have Tre{X^ a ^]) = TrX[ aM ^ = MTrXy a ^y 
More explicitly we introduce the linear real function of domain Dn : RO — ► R, defined 
on any region operator Xa € RO, determined by a domain A e L, of the phase space, 
as the area of that domain. Then we obtain from the last equation the total area of 

M-l 

the M-gon as Dn[e{X [a ^)\ = U #(^r)-[ a > Tt\ = K M L wh ere [a, ff] and [a,M] 

j=0 

parametrize respectively the triangle and the polygon, and the notation R(^0-).[a, jj], 
signifies the area swept out from the clockwise rotation of the triangle [a, by the 
angle ^. 

Let us consider also the dual e* of the CPTI map e, that operates on the density 
matrix and produces the same qpi. We recall that for some p density state operator the 
integral of its associated Wigner function W p (x,p) over the area S a ^M of a a— canonical 
M-gon gives the qpi Q a ,M[ W p\ = Is M w p( x >p)^t- Tlic du& l ma P £ * is defined as 
follows 

QaMW P } - Tr(pX aM ) = (p,X a . M ) - (p,e(X 0iS )) = (e*{p),X a ^) (17) 

The bracket in last equation above stems from the observable-state duality and refers to 
trace class operators, as the density and the region operators are considered to be. This 
last equation determines the dual CPTI map e* , which explicitly reads 

^)=£*(^)V^). (is) 

3=0 

Next we seek to construct operationally the polygon region operator X^^m] , if the triangle 
region operator X^ a ^ is given. Since the maps £*,£, are not easy to be construct 
physically we can use the standard Naimark extension theorem [8] , which applied to CP 
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trace preserving maps allows to implement them by means of some unitary operator in a 
space that is an extension of the original one by an auxiliary (ancilla) space, followed by a 
projection to the initial space by means of a partial trace. This standard unitary dilation 
theorem of a CPTP map, will be adapted here to our case for the CPTI dual maps. To 
this end we introduce an M dimensional ancilla space Ha = {\m) ,m = 0, 1, M — 1} 
and choose some operator |fc) (fc| of it. Then on the total space H <£) Ha, we define an 
operator V, by means of which we write for the maps e, e* the respective expressions 

e{X [am ) = X [aM] = Tr A V{\k) {k\®X [a ^)V\ s*(p) = Tr A V\\k) (k\ ® p)V, (19) 

where we choose (j\ V \k) = R{^f), j = 0, 1, . . . , M - 1. 

In physical terms, this realization describes the coupling of the given quantum density 
matrix p or the quantum observable X^ a ^i, with an ancilla, "atomic" M— dimensional 
system, their dynamical interaction by means of the operator V„ and their subsequent 
decoupling as a result of partially tracing the ancilla system. The physical resources 
of such an operational realization is however very demanding since they require the 
interaction with a high dimensional "atomic" system. More feasible realizations which 
are also motivated geometrically by the various ways of covering a polygon by rotated 
and reflected triangles as desired. In fact one would like to have realization that use only 
two dimensional ancilla systems i.e qubits. Although there is no explicit proof detail 
search shows that there are no realizations that would use only qubits. 

In fact alternative realizations would involve multiple CP maps and Kraus generators [8] . 
Those generators would be the identity map 1, a rotation operator R(<p) of angle ip, or a 
reflection operator A(r) with respect to some axis e.g going through the origin and mak- 
ing an angle r, with the x-axis. Two types of maps operating on quantum observablcs 
are then introduced, one for rotations and one for reflections, which read respectively 

e v {X) = X + R{ V )XR( V )\ e r (X)=X + A(r)XA( r y. (20) 

For observables that are region operators corresponding to a phase space domain with 
parameters specified by D, we denote the region operator by Xjj. Then the linearity of 
CP map yields for the transformed, by rotation and reflection maps, region operators the 
respective relations £ ip {Xr>) — ^Du-R(ip).r>! an d £ r (Xo) = Xu U A( r ).D- The notation D U 
R{<p).D (DUA(r).D), denotes that the domain Bn[e v (X D )}= DL)R(ip).D, (Dn[e r (X D )} = 
DUi(r).fl), of transformed, by the rotation (reflection) CP map, region operator, is the 
union of the domain of the untransformed one by the rotated R((p).D (reflected A(r).D) 
image of it. This geometric property common to all CPTI maps, would probably justify 
the name "copy-paste transformation" for them. Dually we consider next the CP map 
e\ , k = ip, r, with respective Kraus generators S = R(p>), A(r), acting on some density 
operator. The introduction of an ancilla Hilbcrt space Ha = span{\&) , \& ± )}, with 
orthogonal basis vectors i.e ($1$^) = 0, where |$) = a |0)+6 |1) , allows one to implement 
those CP maps unitarily as 



e* k {p) = p + S^pS = Tr A (V k |$) ($| (g) pVt). 



(21) 
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The Kraus generators are obtained by the relations 1 = ($| Vk |$), and 5* = ($" L | Vk |$)- 
The unitary operator Vk, determines given CP map up to a unitary local (operating 
non trivially only in ancilla space ) operator i.e for W a unitary operator in Ha, the 
transformation Vk — ► W <g> 1 Vk, generates the same CP map. In our case such a Vk 
operator in the {|$) , j^" 1 )} basis is 

V k = |$) ($| <8> 1 + |$) ($- L | ® (-5) + |$-L> ($| <g> S*t + O 1 

= E iJ =i,2l* < )<^'l®K i =( <yt if) ( 22 ) 

where {\&) , |$ 2 )} ={|$) , j^)}, and VJ, 11 = ^ 22 = 1, ^ 21 = -F fe 12 = 5. 

Next introducing and combining the next two maps e\, e\, we can construct as an 
example, the hexagon region operator with radial parameter a = This realization 
reads 

s 2 (ei(X [a ^ ] )) = X M , (23) 

where, using the equality R(^-) = R{^Y = R{— f ) (see fig. 2 also), we have the explicit 
maps 

ei{X [am ) = X [am + RC-)X [am RC-y + R(^)X [am R^)\ (24) 
e 2 {X [am ) = X [am + A{^)X [am A{^. (25) 

The domain of the initial isosceles triangle [a, increases by successive "copy-paste" 
induced by the two maps until reaching the hexagon domain [a, 6] as the following equa- 
tions show, 

Dn[ £l (X [am )} = [a, — ] U R(-).[a, — ] U R(--).[a, — ] 

Dn[ £2 ( £l (X [ai¥] ))] = Dn[ e i(X [aiW )]U A(|).Dn[ei(X [0i¥] )] = Dn[X [a , 6] ] = [a, 6]. (26) 

See fig. 2 for the geometric effect of the combined CPTI maps to construct the hexagon 
region operator. An implementation of the maps based on unitary operators is then given 
by the relation 

e 2 (e 1 {X [am )) = Tr A ® Tr A (V r=i v=±f (e 1 ) (e 1 ) ® (D 1 ) (^| 

^ [ai ^]K t =f ¥ ,=±f) 
= E,j=i, a E w =i, a t£L ±f x wm vl k 2 ± , v r % (27) 

where two ancilla Hilbert spaces, one three dimensional i.e H\ = spanjlO 1 } , |0 2 ) , |B 3 )}, 
and one two dimensional H\ — span{|$ 1 ) , |$ 2 }} have been used. The operators in- 
volved are also of two kinds: the first one V v =±il is a three dimensional unitary dila- 
tion of the Kraus generators of the CPTI "rotation" map e\ i.e 1 = (0| V v= ±^ I® 1 ) i 
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i?(§) = (l|n,= ±f le 1 ),^ ) - {2\V v=±i (6 1 ), with the property ^=± f Vj =±f = 
31,. The second one K=§ , is also a two dimensional dilation of the "reflection" CPTI 

map e 2 defined in the canonical basis as V r =* = I \ 2 ) , and is such 



1 

that V r =s.V}_* = 21. It should be noticed that the reflection with respect to the 

2 2 

y-axis i.e (x,y) — > (—x,y), is not a canonical transformation, but since the domain 
Dn[ei(X[ a 22r])], is the semicircle we actually need the reflection with respect to the 
origin i.e (x,y) — > (— cc, — y), which in fact is a canonical transformation, and it is re- 
alized by the parity operator II = e mN = (— 1)^, so we will in effect use the operator 

K=f = ^ 1^ ) ' ^ n< ^ ecc ^ tnc P ar ^y operator with properties lit = n, and II 2 = 1, 

operates on the canonical generators as Ilall = —a, Hatn = —at, and on the displace- 
ment operator as UD a Il = D_ a . 

4. Bounds, Tiling and Majorization for QPMs 

Let us start by pointing out the covariance property of the displaced parity kernel 
operator A(g,p) = D(q,p)HD^ (q,p) of the region operator under the displacement op- 
erator, which reads 

D{q ,p')A(g,p)£%' ,p') = A{q + q ,p + p). (28) 
We now introduce a CPTI displacement map by the equation 

e (X)=X + D(q,p)XDi(q,p). (29) 

Two particular such maps are the q — west and the p — north displacement maps defined 
correspondingly as 

e w , = e , , e N ,=e , . (30) 

9 (9 .0)' P (0,p') v ' 

For a region operator X[ qp ] or -X"[ c ,c;d] w ith support [q,p] a rectangle with sides q,p along 
x,y axes, or [c, c;d] a disk of diameter d and center at the point (c, c), the previously 
introduced domain operator is respectively 

Dn[X [qM ] = [q,p], Dn[X [c . A ] = [c;d]. (31) 

Also we will use the following symbols to denote q —west and the p —north displacements, 
of respective domains of phase space 

W q >.[q,p] =[q'+q,p], N p , . [q, p] = [q, p + p] , (32) 
W d .[c,c;d] =[c + d,c;d], N d .[c, c; d] = [c, c + d; d] (33) 

The CPTI maps that can induce the respective displacement in the domains of rect- 
angle region operators are 

e$(X [qiP] ) = X [qtP] + D(0,p')X lqiP] D\0,p), (34) 

eJ(X [qM )=X M +D(q',0)X M DHq',0), (35) 
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and similarly for the respective displacements of a [c, c; d] disk, 



e^{X [c , c . d] ) = X [c>c;d] +£)(0,p')X [CjCi<fl £)t(0,p'), 
eJ(X [CtC .^) = X [CtC . d] + D(q ,0)X [CtC . d] D\q ,0), 



(36) 
(37) 



Indeed we have for the rectangle that 



Dn[s™,{X [qtP] )} = [q,p] U N p ..[q,p] = [q,p] U [q,p +p], 
Dn[sJ(X [gtP] )} = [q,p] U W q , \q,p] = [q,p] U [q +q,p]. 



(38) 
(39) 



and for the disk 



Dn[e p ,{X [c , c . d] )] = [c,c;d] U N p /.[c,c;d] = [c,c;d] U [c,c + p ;d], 
Dn[eJ(X [c , c . d] )} = [c,c;d]UW q ,.[c,c;d] = [c, c; d] U [c + q, c; d]. 



(40) 
(41) 



In view of fig. 3, we now introduce a rectangle operator X[ a ^], (fig. 3a), by along 
the x axis to west, (fig. 3b), and one displacement of the resulting rectangle by e„ along 
the y axis to north, (fig. 3c). The combined displacements define the map £( a ,/3) realizing a 
single step in the tiling by operating on the square and disk region operators respectively 
as follows 



where the notation X[ c+d . c+d . d ] means the cluster of disks of diameter d, with center at 
the points {(c, c); (c, c + d) \ (c + d, d); (c + d, c + d)} . This notation extended in the form 
X[ c +md.c+md:d], where m — 0,1,2,... , describes the cluster of all disks of diameter d, 
"under" the disk with center at the point (c + md, c + md). 

Referring to fig. 3 (fig. 4) we introduce a west-north tiling of an initial rectangle (disk) 
region operator X^ a ^, (X[ CiC;c q), by successively operating on it with the CPTI maps 
{£( Q ,/3)>£(2 a ,2/3),£(4a,4/3), ■•■}> iU(d,d) , £(2d,2d) . £{id,id) ' -})' wnicn results into the sequence 
of region operators {^[ a ,/3]>^[2a,2/3],^[3a,3/3]i ■■■}, {{X[ c . c - d ], X[ c+dc+d . d ], X[ c+2d . c+ 2d-.d], ■■■}), 
the sequence of domains of which corresponds to a rectangle (disk) shaped west-north 
tiling of the upper left part of the plane. In other words we have the iteration 



for the two types of tiling of the original rectangle and disk region operators. 

Before proceeding let's us mention that the CPTI maps can be unitarily extended 
by utilizing an auxiliary two dimensional Hilbert space and a partial tracing as follows 
(y ^rectangle, disk), 



e^(X [y] )=Tr A (Uj\0)(0\^X [y] U^), e${X [v] ) = Tr A (U^\Q)(Q\ ® X [y] U^). (46) 



X [2a,20] — £(a,f3)(X[ a Jj]) = E fi OE a (X[ a fj ]), 

X[ c +d,c+d-d] = £(d,d)(X[ C:C . d }) = e d o (X[ CtC . d j), 



(42) 
(43) 



[c+md,c+md:d] = e ((m- l)d, (m-l)d) (^[c+(m-l)d,c+(m-l)d;d]) > 



me {2,4, 6,...}, (44) 
me {1,2,3,...}, (45) 
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The unormalized unitary operators chosen for the extension of the west and north dis- 
placement maps are respectively the following two 

T jw_( 1 -D\q,Q)\ N _( 1 -D\Q,p)\ 

«' ~\D(q',0) 1 U p'-\D(0,p) 1 J' ^> 

Let us turn now to the question of determining the quasi probability masses (qpm) 
identified with the eigenvalues of rectangle (disk) region operator X^^^] 
(X[ c+mc i .c+md-d]) at some time [2/i,2v] (m), of the tiling iteration process, from its previ- 
ous one [fi, v\ (m — 1). This would be an important relation for determining qpm's, and 
their bounds i.e the maximal and minimal ones among them, in the course of rectangle 
(disk) tiling, once the qpm's of the initial rectangle (disk) region operator is given or 
can be measured or estimated. To this end let us take the canonical decomposition of 
the matrices of a region operator X = V'X d V, before and after a general tiling CPTI 

map i.e of the operator X = e , , (X) = W^x"'w. Here V, W and X d , X*' , are the 

(i ,p ) 

respective diagonalizing unitary matrices and the diagonal matrices of the eigenvalues 
appearing in the canonical decomposition. The latter two are related by the equation 

X dl = W[V^X d V + D(q,p)V^X d VD\q',p')}W^ (48) 

Let A' = = ({X d ')jj)°^ 1 , the eigenvalue column vector, and similarly from the 

eigenvalues of the X d the vector A. In view of the last equation these two real column 
vectors are connected by the following matrix E 

A' = (WV^ o Wt + WD(q , p ')V ] o WD(q',p')V^ A = SA. (49) 

Above, the overbar denotes complex conjugation and, the elementwise or Hadamard 
product between matrices defined as (AoB)ij = AijBij, has been used. Since for a given 
unitary matrix U the Hadamard product U o JJ is a column and row stochastic matrix 
(bistochastic), i.e for the unit column vector e = (1, 1, ...) T , we have that U o Ue = e, 
and e T U o U = e T , then our present sigma matrix has column and row sums equal to 
two i.e Se = 2e, and e T S = 2e T . 

In order to apply these facts to two successive region operators in the rectangle tiling, 
we introduce first their canonical decompositions 

*[2„,2„] = W^Xf^W = WUiag{\ [2llM )W, (50) 
X M - V^Xf^V = VUiag(X M )V. (51) 

Then we obtain the following relation among their diagonal components 

X[2vM = W{V^Xl M V + D(^Q)V^Xl M VD(^ 

+D(0,v)ViX^ v] VD(0,v)i 
+D(0, v)D{^ 0)V*X^ v] VD(n, 0)t£>(0, ^}W^ (52) 
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The diagonal elements of these matrices are the column vectors A[ M) „], Ap^i/], and they 
are identified with qpm's before and after the tiling map. They are explicitly related as 
follows 

\2„,2u] = (WV^ o Wt + WD(p, 0)0 o WD(n,0)V 

+WD(0, v)V ] o WD(0,v)Vl + WDfa, v)V^ o WD{fi,v)V^ X M 

= TX[^]. (53) 

The matrix gamma introduced above updates the vector of qpm's in the iteration of 
rectangle tiling, and it has column and row sum equal to four i.c Te = 4e, e T T = 4e T ,. 
This last fact signifies the quadrupling of a initial rectangle at each tiling step. As 
to the disk tiling, the analogous relation results after simplifying by setting W = 1, 
since the disk range operator X[ cc . d ] is diagonal in the number state basis i.e X[ cc . d ] = 
X[cc-d] = diag{\[ c ^ c -d]) ■ Then the relation connecting the column vectors \ c ,c;d] an d 
\c+d,c+d;d]i of the diagonal elements of the respective disk range operators -X"[ c +d,c+d;d] = 
£(d,d){ x [c,c;d]) = VUiag(X [c+dtC+d . d] )V, and X [c ^ c . d] , i.e the operators before and after a 
disk tiling map, reads 

\c+d,c + d-4] = (V ° + D(d, 0)0 o D(d,0)Vi 
+D(0, d)0 o D(0,d)Vi + D(d, d)0 o D(d,d)V^ A [c , c;d] = EX [c . c . d] . (54) 

As with the matrix T, the matrix E has column and row sums equal to four. 

Remarks: 1) The eigenvalues Ap^i/] i-e the qpm's, will result if we evaluate the 
matrix elements of some square region operators in some of their eigenvectors. These 
eigenvectors however may not be known and not even been relevant vectors for some 
problems, in such cases the above recurrent relations may not be useful. Instead the 
determination of matrix elements in the number states which are physically relevant 
states, is desirable and this will be pursued next. 

2) The same remark is in fact true for the determination of eigenvalues X[ c+md c+md . d ^ 
i.e the qpm's in the case of disk tiling. So in the following we evaluate the diagonal matrix 
elements of the respective region operators in the number states. 

Lemma 4. Let the transformation D — > e(D) = X = J A V(x)DV^ (x)dfi(x), of the 
operator D = X)fe=o l dk\k)(k\, diagonal in the number state basis, and V = V(x) a 
unitary operator depending on some variable x with values in some domain A, on which a 
integration measure d/j,(x) is defined. Let the vector \^) = J2k=o l dk\k), then the diago- 
nal elements of operator X, are determined by (k\X\k) = (k\ ^ J xeA V(x) o V(x)dfi(x)j \^) . 

The proof is straightforward, and the lemma is also true for the case of transforma- 
tions defined by discrete sums. All the above can now be summarized in the following 
proposition: 
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Proposition 5. Let the rectangle and disk tiling with iteration maps £( M) „) and 
£ (md.md) respectively, acting on rectangle and disk region operator as e^^(X^^) = 
X[2n,2v] an d £(md,md){ x [c+md,c+md-d}) = X [c+(m+i)d.c+(m+i)d-d] ■ Let the canonical de- 
composition of the rectangle tiling range operator X^,^ = V^X^ = V' f diag(X^ . V \)V, 
with \n >v \ = {^{^yjIfeLo' ^ e c °l um n vector of its eigenvalues. Let -X"[ c + TO d,c+md] = 
V^X* c+md c+md ]V = diag(\[ c+md . c+md ])V been also the canonical decomposition of 

the disk tiling range operator with \ c +md.c+md] = {^[ c +md c+md]}fe=0 ^ e column vector 
of its eigenvalues. Given that the eigenvalue-vector problem has been solved for these 
operators after some given tiling step, then their diagonal elements in the number state 
basis, (which physically signify the quasiprobability mass of the Wigner function in the 
respective number states), are given by the two corresponding expressions, 

(fc|-X"[2„,2,#> = (k\ \v*oV + D(jt, 0)0 o VDijt, 0)t + £>(0, o VD(0, i/)t 



+D(n,v)V* o VD{^\ (55) 

and 



{k\X [c+[m+1)d , c+(m+1)d . d] \k) = (k\ oV + D(md,0)V^ oVD(md,0y 



+D{0,md)V^ o VD(0,md)i 



+D(md, md)V^ o VD(d, md)tj \^ [c+md . c+md] ) . (56) 

where \$ M ) = E fe =o,i,.. A [k#> and ' [c+md,c+md]) = E fe =o,i,... X [c+ md ,c+ m d]\ k ) ■ 

The description of the tiling process that has been given above with rectangle and 
disc tiles is based respectively on the one step CPTI maps £(^ tV ){X^^) = X[ 2fl .2v], and 
£(md.,nd){X[ c+m d.c+ m d-d\) = x [c+(m+i)d,c+(m+i)d;d]- These maps induce in the eigenvalue 
vectors of the region operators the corresponding transformations X[ c +d.c+d-,d] = EX[ c ,c;d], 
and X[2n,2v] = TA^ ^j . The matrices E and T involved, are almost bistochastic, since as 
we have shown they have column and raw sums equal to four. This property allows to 
estimate the behavior of upper and lower bounds of the eigenvalues of the region operators 
after a one step operation of the tiling process. Specifically in the subsequent corollary we 
will determine the relation before and after the performance of one step tiling, between the 
extremal values of the eigenvalues of rectangle and disk region operators. These values 
of course stand for the extremal values of the quasiprobability mass of the respective 
region operators in their associated eigenstates. These extremal values provide useful 
information by means of which we can estimate the interval of values of the qpm's in the 
course of tiling process. 

To this end let us introduce the following notation: let the matrix S which stands 
for the matrices E and V used above, sharing the same properties. This matrix will be 
used to write generically A = SX, for a transformation of the eigenvalues in a one-step 
tiling. Let us introduce the non- increasing (non- decreasing) ordering of the eigenvalue 
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vectors i.e respectively A^ = (X[, \$, , A| +1 , ...), with Xj > Aj +1 , i — 1,2,... , and 
AT = (Al,A|,...,AT,AT +1 ,...), with A] < Aj +1 , i = l,2,... . 

Corollary 6. Let the transformation A = SX, of the column vector of rectan- 
gle/disk region operator eigenvalues in a one-step tiling map. For the non-increasing 
and the non- decreasing eigenvalue vectors the following respective majorization relations 
are valid, 4A^ >- X^' and A^/ >~ 4A^. Also for the maximum and minimum values of 
the eigenvectors before and after the transformation the "squeezing" property is valid 
4A m ; n < A min < A max < 4A max . 

' l > i 

Proof: We start with the equation A = SX, but it suffices to assume A — SX , and 

A 1 = SX 1 , for we could always consider two suitable permutations P, Q such that A = 

A 1, i 

PA , A = PA. The equation for the non-increasing ordered vectors A = PSW X , 
involves the matrix PSW T with the same property as the S, matrix i.e PSW T e — PSe 
= 4Pe = 4e and e T PSW T = e T SW T = Ae T W T = 4e T . Similar result are obtained 
for the non-decreasingly ordered vectors. We start now with the equation A = SX 

A/ c I I 

and write A = 4-|A = 4PA , where H is a bistochastic matrix. The transformation 

A^ = HX l , leads to the majorization A^ >- X 1 [10]. Then the relation ^A = X 1 , 

leads to the majorization 4A^ y X^' . The proof of the second majorization relation is 

similar. Finally, taking the maximum and minimum values of qpm's before i.e A max = 

it it 
A 1; A m in = A{, and after i.e A max = A 1; A m ; n = A{, of the tiling map, and using the 

inequalities ensuing from the majorization relations we obtain the "squeezing" property 

4A m in < A min < A max < 4A max B 

5. Conclusions 

Region operators (RO) associated with integrals of Wigner quasiprobability function 
over domains of phase space provide an example of new type of quantum mechanical 
observables. These operators can be considered as a quantum version of characteris- 
tic functions of classical regions, or as generalized operator valued probability measures 
(OVM). The theory connecting geometric transformations of regions to spectral charac- 
teristics of ROs, has been initiated in this work. Applications of ROs and their associated 
qpm's with their respected bounds, to the field of tomographic reconstruction of Wigner 
function over certain regions [2, 3], should be a testing ground for the theory. Also the 
application of region operators in the study of coupled bipartite quantum systems, could 
provide a framework for investigation of the geometric manifestation and quantification 
of quantum entanglement, cast in the language of generalized operator valued measures 
on phase space. We aim to return to these matters elsewhere. 
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Figure captions. 

Fig.l Point and line region operators K = U and Kq = (II) respectively. 

Fig. 2 Construction of hexagon region operators -X"[ a ,6] (2c), by implementation of two 
CPTI maps on the isosceles triangle region operator X[ a * j (2a). 

Fig. 3 Construction of west-north tiling by successive application of CPTI maps on 
the rectangle region operators ^r Q ,/3] (3a). 

Fig. 4 Construction (i) of west-north (4b)and (ii) of reflection-rotation (4d, 4e) tiling 
by successive application of CPTI maps on region operators of clusters of disks (4a) and 
rectangles (4c) respectively. 



